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The 6j-symbol is a fundamental object from the re-coupling theory of SU(2) representations. In 
the limit of large angular momenta, its asymptotics is known to be described by the geometry of 
a tetrahedron with quantized lengths. This article presents a new recursion formula for the square 
of the 6j-symbol. In the asymptotic regime, the new recursion is shown to characterize the closure 
' of the relevant tetrahedron. Since the 6j-symbol is the basic building block of the Ponzano-Regge 

^SJ , model for pure three-dimensional quantum gravity, we also discuss how to generalize the method to 

►..^ ■ derive more general recursion relations on the full amplitudes. 
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The theory of SU(2) representations and re-coupling is known for quite a long time. The main motivation was the 
quantization of angular momenta, with obvious applications in spectroscopy, atomic/molecular physics. Nevertheless, 

^ the study of SU(2) re-coupling coefficients remains an active field of research in modern physics One reason is 

the appearance of a new notion to describe some phases of matter: topological order Its emergence has been 
~T described as a condensation of spin networks ^] (possibly with a quantum group deformation). It also gives a new 
way to think of fault-tolerant quantum computations 
I ■ In the same time, a similar formalism has led to interests in re-coupling coefficients to capture aspects of quantum 
geometry and possibly quantum gravity. The main idea in that approach is that there is no fixed background geometry, 
but geometry is instead a fully dynamical object. Spin foam models have emerged as direct applications of that idea 
(see e.g. Q for a review). A useful model to explore the spin foam program is the Ponzano-Regge model [1, 0] 
which aims at defining the path integral for three-dimensional pure gravity. The latter is a topological field theory: 
there are no local degrees of freedom due to a large set of specific gauge symmetries. Thus some interests in that 
model also comes from topological invariants, and especially quantum invariants since a quantum deformation of the 
' Ponzano-Regge model gives the well-defined Turaev-Viro model Q . 
^T) . These models are constructed by triangulating the space-time manifold and defining the corresponding amplitude 
' as a product of local amplitudes attached to each tetrahedron of the triangulation. The basic building block of the 
] Ponzano-Regge model associated to each tetrahedron is the Wigner 6j-symbol of the SU(2) re-coupling theory. 
• The 6j-symbol is certainly the most studied Wigner coefficient, and we refer to for modern reviews describing 

T-H ' numerous aspects and to for textbook results. An important regime of such a quantum object is the semi- 
^ , classical limit, when angular momenta become large. The asymptotics of the 6j-symbol has been studied and proved 
. ^ ' in [T2| - [l5t , and next-to-leading orders have been obtained in [343 ■ The asymptotics of the Ponzano-Regge model 
, on handlebodies (with more than one tetrahedron) is presented in [l^. Recent works have also focused on more 
■ complicated objects, like the 15j-symbol and EPR amplitude in [20j. An extended Born-Oppenheimer approximation 
- - - has been developped and applied to the asymptotics of the 9j, 12j and 15j-symbols with some large and small angular 
momenta [2l[. One of the author has also studied this regime for arbitrary 3nj-symbols in with a simpler method. 

An efficient way to study such asymptotics as well as algebraic properties of re-coupling coefficients is the use of 
recursion relations on the angular momenta. It is actually a way to define the 6j-symbol [Tl|, [l^ , and some insights 
into asymptotics and recursions on arbitrary 3nj-symbols are given in [23[ . From the point of view of topological field 
theory, recursion relations have been understood as a way to encode (some features of) the dynamics, i.e. characterizing 
the invariance under gauge symmetries [l^, . AH these considerations actually also apply to the four-dimensional 
case f26| and spin foam models for 4d quantum gravity. In this case, we do not yet have the symmetries and the 
dynamics under full control and we strongly think that the development of recursion relations for the spin foam 
amplitudes will be a very useful tool to understand further the structure and invariance of 4d spin foams. 

There exist several ways to derive recursion relations for the 6j-symbol. The original method is to start from the 
Biedenharn-EUiott (BE) pentagon identity. The latter encodes at the algebraic level a topological invariance of the 
Ponzano-Regge model under a special move. It is an equality between the product of two 6j-symbols (the weight for 
two tetrahedra glued along a triangle) and a sum of products of three 6j-symbols (which correspond to three tetrahedra 
sharing an edgejj. Recursion relations are obtained by specializing the BE identity with well-chosen boundary angular 
momenta Il3l. 
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^ This states the invariance under the 3-2 Pachner move. Topological invariance also requires invariance under the 4-1 move which is 
actually induced from the Biedenharn-EUiott identity together with the orthogonality relation of 6j-symbols. But it is only formal due 
to divergences [27l. |2S(I . 
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Another path towards recursion relations is to think of Wigner symbols as evaluations of spin network states. 
Then, one can interpret recursion relations as Wheeler-DeWitt equations in the spin network basis, generated by 
some Hamiltonian operators which annihilate the evaluations. Such an operator is expected to be the generator of a 
gauge symmetry [24-:2|jj . That is a way to see how recursion relations encode and provide the spin network states of 
Loop quantum gravity with their dynamics. The standard recursion on the 6j-symbol was derived in that way in f25| , 
where the Hamiltonian operators form a first-class algebra classically, which means that they really generate a gauge 
symmetry. 

The final method which we would like to report on here is to start with the expression of the 6j-symbol as an 
integral over four copies of SU(2) (see e.g. That group integral approach was already proposed in an earlier 

study 24]. In this previous work, the method could only be applied to a special class of isosceles 6j-symbols (and 
to the lOj-symbol of the Barrett-Crane model for 4d Euclidean quantum gravity). The idea is that measures in such 
group integrals have some specific properties which can be used to get recursion formulae. Further, those measures 
have a geometric nature which eases the interpretation in the asymptotics. 

However, the case of generic 6j-symbols has a more complicated measure. Still, it is possible to extract recursion 
formulae using similar techniques, as we show in the present paper. The organization is the following. In the Section 
m we present the integral which is our starting point. In the Section|lTl we prove a new recursion, acting on the square 
of the 6j-symbol, whose asymptotics and semi-classical interpretation are discussed in the Section iHll Finally we open 
in the Section IIVI the interesting possibility of deriving recursion relations in a systematic fashion for arbitrary spin 
foam amplitudes in the Ponzano-Regge model, as opposed to a single tetrahedron. 



I. 6J-SYMBOL AS A GROUP INTEGRAL 



Let us consider a tetrahedron with its six edges labeled by i = 1, . . . , 6 and its four vertices labeled by a = 1, . . . , 4. 
We attach an irreducible representation of SU(2) to each edge, i.e we associate a spin ji G N/2 to the edge i. 
Equivalently, we will refer to these spins as j(ab) where the (symmetric) pair of vertices a, b uniquely defines the 
edge i. At each vertex of the tetrahedron graph, we associate the unique (up to normalization) 3-valent intertwiner 
between the representations living on the three edges attached to that vertex. The intertwiner is the invariant vector 
in the tensor product of the three representations meeting at the vertex. Intertwiners are given explicitly by the 
Clebsch-Gordan coefficients, or more precisely by the Wigner 3jm-symbols. Finally, tensoring these four intertwiners 
living on the four vertices of the graph and taking the trace on magnetic indices, we obtain the 6j-symbol. 

The square of the 6j-symbol can be expressed as an integral over the group SU(2), using the integral expression of 
the products of 3jm-symbols 

i—l ^ ^ ^ ^ 

where I?(^'^(g) is the Wigner matrix of 5 G SU(2) in the representation of spin ji and the quantities into brackets are 
the 3jm-symbols. This way, we get one integral over the group for each vertex of the graph. Taking the trace leads 
to an expression for the square of the 6j-symbol as an integral over SU(2)"' (see e.g. [l4| for more details), 

where the integral is taken with the Haar measure \dg\ over SU(2) and the character xAd) = Tr^Z?'^-') (5) is the trace 
of the Wigner matrix 13^^' (5). This is the starting point of our analysis. 

Since characters are invariant under conjugation by SU(2) group elements, we can write the above expression as an 
integral over the class angles Qah of the group elements gag^^ ■ The change of measure was computed in p^ . 

= ^ / [ddab? \ ^ . n «in(2ja6 + l)eab , (3) 

TT Jd ^det4x4(cOs6'afc) (^^j 

where we take the convention that 0aa = and thus cos 9 aa = 1- The measure is given in term of the determinant 
of a 4 X 4 Gram matrix (cos^ab). It has a simple geometrical interpretation: •\/det(cos dab) gives the volume of the 
(dual) tetrahedron on the sphere S^, whose spherical edge lengths are the angles dab- The above formula can be used 
to extract the asymptotic behavior of the 6j-symbol at large spin through a saddle point analysis JjJ . 
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The domain of integration D is the set of aU possible spherical tetrahedra, i.e. the subset of [0,7r]^ satisfying 

dac < dab + Obc and Oab + Obc + Oca < 27r. 

Since it will be our departure point, let us sketch the derivation of ([3]). Among the four group variables in ([2]), 
we can absorb, say, (74 in a re-definition of the others, ga QaQ^^ for a = 1,2, 3. Clearly, that does not change the 
arguments of the characters. Moreover, due to the right invariance of the Haar measure, dgi, dg2, dg^ are unchanged, 
and the integral over 34 becomes trivial. Then we write each group element ga — exp^iOaifia ■ (j), where Oa^ is the 
class angle, and fia £ the axis of the rotation, and the measure reads 

2 

dga = — Sin^ OaidOaA d^ha ■ (4) 
IT 

The class angles are kept as integration variables since they are the arguments of the characters with spins ja4- As 
for the six variables {na)a=l.2,3^ we split them into two sets. One consists in the dot products Uab = fia ■ fib, and the 
other set form a global rotation (three components) which leaves the variables Uab invariant and trivially factorizes. 
Finally, we notice that the arguments of the characters with spins (jab) with a, 5 7^ 4 are class angles, which we denote 
9ab, which are related to the variables (uafc) by 

cos dab = cos da4 COS ^{,4 + sin 9a4 sin 9b4 Uab ■ (5) 

The final step is to change the variables from Uab to 6ab- The Jacobian is straightforward to evaluate and leads to 



II. NEW RECURSION 



Let us introduce the following notations, 

P[Xab]=det{Xab), and T^Jijab) = fijab ± ^ 



(6) 



P is the determinant of a collection of numbers (Xab), and is polynomial of degree four in the six variables Xab = Xba, 
Xaa — 1- The operators shift the spin jab for a given pair (ab). The new recursion we are going to present is 



'iP[-2iT^b + Tab)] ((2jo^ + 2)r+, ~ 2j„^r-^)+9^„,P[i(r++Tj](T+ -T-^)2] {6jp = 



(7) 



for each symmetric pair (a/3) of vertices, with a, /3 = 1, 2, 3, 4 and a ^ 13. 

Our method is fairly simple and works exactly like Schwinger-Dyson equations in field theory. We multiply the 
integrand of ([3]) with a well chosen insertion and take a total derivative and integrate over all angles. 

Singling out one edge (a/3), we consider the following integral 



[dOab]^ de^p sin 6*0/5 . /det(cos 9ab) J| sm{2 jab + l)9ab 

(ab) 



0, 



(8) 



which is indeed zero since the factor det(cos6'a6) vanishes on the boundary of the domain D 'T^. 

To expand this integral, we use the notation — det4x4 Xab- The measure term is given by its evaluation on 

Xab = cos 9ab- Wc computc explicitly the derivative in the previous expression and distinguish the two terms where 
we differentiate the Gram matrix determinant or the other terms. 



= J [d9ab] 



det(cos Oab) 



^det(cos 9ab) 
9x^^P[Xab]\x=cose .2 



(cOs9ap Sm{2jap + l)0al3 + {2jal3 + 1) siu cos(2j„/5 + l)0al3 



2^det(cos6' 

ab 

which can be re-written as 

det(cos 9ab) 



sin 9a/3 sm{2jai3 + 1)6*0/3 



Y[ sm{2 jab + l)9ab 
J (a6)#(a/3) 



(9) 



= / [d9ab] 



Y/det(cos 9ab) 

dx^,P[Xab]\x=, 



+ 2) sin(2jQ^ -I- 2)9ap - 2jafi sm2 japOap 



v/det(cos 9ab) 



(1 - cos 9al3) Sm{2jal3 + 1)0q/} 



Yl sin(2j„b + l)9ab ■ 

{ab)^(afi) 



(10) 
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Now all the cos 6'- factors can be re-absorbed as shifts in the spin labels jab- Indeed, the multiplication by a factor 
cos 9ab corresponds to the operator ^(r^ + T^) due to the trigonometric identity 



cos fab Sm{2 jab + l)Oab 



Sm{2jab + 2)dab + sm{2jab 



^ab + Tab gi^(2j,;, + l)6ab ■ 



2 2 
That leads to the following recursion relation on the squared 6j-symbol 



(11) 



{6j}2 = 



. This is ([7]) . The non-triviality of the relation is ensured by the presence 
of the factors (2j + 2) and 2j in the first line, which implies that the coefficients of this equation can not all vanish. 



where we used that T nT „ 

ajj ajj 



To better understand the explicit structure of this relation, we need to expand the Gram determinant, 
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P[Xab 



E 



(12) 



where we sum over all permutations of 4 elements. We classify them by their number of fixed points. The only 
permutation that fixes all elements is the identity. Then we get the two-point cycles, which exchange two elements 
without touching the other two elements. These contribute as —X^^^ to the determinant, where a and b are the 
two elements which are interchanged. Then permutations with a single fixed point are given by the 3-cycles, e.g. 
1 — !■ 2 — ^ 4 — > 1 leaving 3 invariant. These give terms of order 3. Finally, we have the six 4-cycles, which lead to term 
of order 4 in the polynomial. More explicitly, we can compute 

-I-2X23X24X34 + 2X13X34X14 + 2X12X24X14 -|- 2X12X23X13 



I v2 "v^2 I v2 "v^2 I v2 -y 
+^12-^34 + ^13-^24 + -^14^: 



23 



(13) 



—2X12X23X34X14 — 2X13X34X24X12 — 2X14X24X23X13 , 
with 4! = 24 terms. We also compute its derivative with respect to, say, X12, 

dxiiP — 2 (Xi2(X|4 — 1) + X23X13 + X24X14 — X23X14X34 — X24X13X34) . (14) 
In the end, we obtain a recursion relation which involves mixed shifts on the spins jab up to five ±i-shifts. 



III. ASYMPTOTICS AND SEMI-CLASSICAL GEOMETRY 



A. The Recursion as a Closure Condition 

The most well-known (and useful) recursion relation on the 6j-symbol is a second order difference equation [ll[, 

A+i[j.] + ^ + Mj.] + ^ ^ = 0, (15) 

^ ' \ J4. J5 jej ^ ' \j4 J5 Jej ^ ' \ J4 J5 JeJ ^ ^ 

and there also exists a similar second order recursion with ±i-shifts. The explicit coefficients can be found in [ill [Tsj 
for instance. This equation can be used to derive all 6j-symbols from an initial condition, and is thus used numerically 
to generate lists of values of the symbol. This means that in principle the recursion relation ([7]) could be obtained 
from this one. This seems however to be quite complicated. 
Still it is interesting to compare both: 

• the standard relation only has three terms, but complicated coefficients, 

• while our new recursion relation has many more terms (corresponding to the various combinations of shifts), 
but all coefficients are trivial (either ±1 or ±{2j + 2), or ±2j). 
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Hence it certainly can not be used in efficient numerical simulations, but in contrast its geometrical interpretation at 
leading order is straightforward, as we will now explain. 

In geometric terms, the recursion (llSp has been interpreted in [2^ as the statement of invariance under an elementary 
deformation on a tetrahedron whose result is a small displacement of a vertex. This simply comes out of interpreting 
some specialization of the Biedenharn-Elliott pentagon identity [ll| in the 3-2 Pachner move. From this point of 
view, the recursion ([7]) should certainly be seen as a non-trivial combination of such elementary moves, successively 
performed. The question is then what geometric properties are encoded into such combination of moves. While the 
standard recursion (|15p encodes the full symmetry of the model (topological invariance), the new recursion should 
only probe some simpler geometric properties of the 6j-symbol. To extract this information, it is far easier and more 
fruitful to look at its asymptotic form. At the leading order, this is just the closure of the corresponding tetrahedron. 

The asymptotics of the 6j-symbol when all spins are homogeneously scaled i^ (see [gI [i^ - [T^ ) 

Here, V[ji\ is the volume of the tetrahedron with edge lengths ji + ^ ^^^"^ = J^iUi + 5)</'jbi] the Regge action 

for the same tetrahedron. The angles arc the dihedral angles of the tetrahedron computed as functions of the 

edge lengths. In particular, they satisfy the following property, 

det(cos(/.a6) = . (17) 

4x4 

Mathematically, it means exactly that the angles (f)ab are dihedral angles for some tetrahedron. Equivalently, we will 
refer to that condition as the closure of a tetrahedron. This is part of the information on the 4>ab[ji]- The other 
information that is needed to recover the asymptotics is that they are the dihedral angles for the tetrahedron whose 
lengths are {jab + 1)- Those two conditions are the roots of the first order formulation of Regge calculus, where lengths 
and angles are independent variables [1^. They can be both extracted from the recursion ()15p . fl3l [. 
Squaring the asymptotic approximation, we easily get 

{6j}' - 2A^V[j{\ ' ''''^^ ^""'^''^^^ ' ^^^f^'^^ " ^^^^'^ ^^'^'f^'^^ ■ ^^^^ 

This leading order asymptotics goes as as we scale the spins homogeneously. 

Let us insert it in the new recursion relation and check that it solves it at leading order in the large spin limit. 
Ignoring the volume pre-factor (which only contributes to higher orders), we have two terms in the asymptotics: the 
constant term 1 and the oscillatory term sin{2 Sji[ji\). 

The constant term trivially solves the recursion relation because of the perfect balance of plus and minus signs 
in the equation. More precisely, both the shifts due to the Gram determinant (since e(cr) = 0) and due to the 
(1 — cos^ 6) factor (and actually also the shifts due to dxP) vanish. 

The non-trivial check thus comes from the oscillatory term. Let us look at the behavior of sin 2Sr when spins are 
shifted. By applying basic trigonometry, we easily see that 

{T^b+T-,) sin(25fl[j-]) = 2cos0,6 sin(25«[j,]) , 

- T-J sin(25fl[j-]) = 2 sin(pab cos(25fl[j,]) . (19) 

Keeping only the leading contribution to the recursion relation, i.e. only the terms with an explicit jap, we see that 
it holds if and only if 

P\\{T^b + T-,)] {T+p - T-p) sin{2SRm ^ sin0„^ det cos</>a6 cos(2^^[j,]) = . (20) 

This is obviously true due to the closure of the tetrahedron, which directly implies the vanishing of the Gram 
determinant of the geometrical dihedral angles, det4x4 cos 0af) — 0. 

Equivalently, if we try to find solutions to the new recursion a la WKB, we find that for any oscillations going like 
exp(z ^ Jab '^af)[jj]), the angles uJab are constrained to be dihedral angles for some tetrahedron, whatever they are as 
functions of the spins. 



^ Note that it can be efficiently derived from |(7]l. 
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If we want to push the analysis to higher orders, we have to include the next-to-leading order the 6j-symbol 
(e.g. [l^ and the corrections due to the volume pre- factor. The geometrical interpretation of the higher order 
corrections to the 6j-symbol itself is not clear at all, so it might be enlightening to push this analysis and see if our 
new recursion relation can provide a natural geometrical meaning. The hope is that our recursion does not encode 
the full symmetries of the 6j-symbol but focuses on its closure which maybe partially characterize the next-to-leading 
orders of the 6j-symbol. 



B. Comparison with other Closure Equations 



1. Asymptotic Equations 

One of the key idea in loop quantum gravity and spin foam models is to make sense of quantum geometry using 
SU(2) re-coupling. From this view point, it is expected that the relevant SU(2) objects satisfy such closure equations. 
So far, it has been mainly investigated in the context of models for four-dimensional quantum geometry. It has been 
recently shown in [26| that closure relations generically holds asymptotically in the large spin limit when the leading 
order oscillates with the Regge action of the 4-simplex. 

Consider a c?-simplex (with c? = 3,4) with its {d — 2)-subsimplices labeled by spins (jab)- Interpreting the latter as 
the volumes of the (cZ — 2)-simplices, we assume we are in a regime where a classical geometry of the simplex with those 
volumes is available. This assumption is obviously satisfied for the 6j-symbol, but is more subtle in four dimensions 
since ten areas do not always determine a flat 4-simplex. It can be reached nevertheless by using an overcomplete 
basis of intcrtwiners which typically provides the normals to the triangles in addition to their areas. The assumption 
then implies that dihedral angles are well-defined and depend on the spins. They may also depend on the normals 
but we can neglect their variations to leading order. 

If the amplitude goes like Vl^[jab] ^ A exp{iSR[jab]), then it satisfies the following equation for each spin jab, 

{{T^bf + (Tab)') Wijab] - 2 cos Cj^ab W[jab] ■ (21) 

It means that multiplication by the cosine of a dihedral angle can be compensated by shifts on the corresponding 
spin of ±1. Then taking the antisymmetrized sum over all permutations of products of that equation, we form on the 
right hand side an identically vanishing factor, det (cos (j)ab) — 0, so that 

P[{T^,f + iT-,f] W[jab] = . (22) 
That provides a definition of a closure operator. 



2. Exact Equations 

The above relation only holds asymptotically, for the 6j-symbol, the 15j-symbol and more complicated objects f2^. 
Hence getting an exact closure relation for the 6j-symbol is an interesting achievement. 

The technical difference between and the asymptotic form of the exact equation ([7]) is the following. The 
operator (TJ^)^ + (r^^^)^ which is the argument of P in (1^^ generates shifts of ±1 while in ([7]), the argument of P 

K^oh + ^af)) which generates shifts of ±i. That difference is compensated since P acts on the 6j-symbol itself in 
but on its square in the exact recursion. Note also that the exact recursion picks up additional terms in the 
sub- leading orders, due to the derivatives of P in 

The new recursion ([7]) is actually the second known example of an exact relation whose leading asymptotic behavior 
reduces to the closure of a simplex. The first example we have found is the lOj-symbol. The similarities of the derivation 
with the present situation strengthens the potential generality of the method. Hence we briefly sketch a few steps of 



that derivation which was reported with details in 2J| . 



The lOj-symbol is a quantum weight associated to a 4-simplex in a model for four-dimensional geometry known as 
the Barrett-Crane spin foam model (soj . It is a function of ten spins j(ab) living on the 4-simplex boundary graph 
(the complete graph with 5 vertices, where each vertex is linked to the other four vertices by a single edge) and its 
asymptotics displays oscillations with the Regge action of the 4-simplex [31]. It is simply expressed as an integral 
over five copies of SU(2), 



(23) 
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Once again, we can write this as an integral over the class angles of the group elements gaS^ and the structure of 
the lOj-symbol is actually simpler than that of the squared 6j-symbol [l^ : 



{lOj} = 




where we simply have to impose that the 5x5 Gram matrix has a vanishing determinant. This constraint states that 
the ten angles Oab are the dihedral angles of a true geometric, flat 4-simplex. 

To get a recursion, we insert the determinant of the Gram matrix in the integral, which then vanishes due to the 
constraint in the measure term, 

/ [dOabY^ 5{det{cos 6 ab)) det(cos0„b) ]T sin(2j,b + l)6ab - . (25) 

/ 5x5 5x5 ^ 

{ab} 

The new factor det5x5(cos ^afj) is a polynomial in the variables cos^af,. Each cosine actually produces a ±i-shift in 
the ? o>i's spins labeling the lOj-symbol, due to the trigonometric identity (ITT|) . 

In [2J], it was checked that amplitudes going like exp{iJ2(ab)('^ja.b + ^)4'ab), i-e. which oscillates with the Regge 
action of the 4-simplex whose triangle areas have the value^j {jab), do indeed satisfy the above equation. That is 
because, not very surprisingly, it reduces to det(cos i/iab) in the asymptotics. 

The surprising fact about that equation on the lOj-symbol is that it is exact although the Regge oscillations are 
definitely not the dominant contributions to the asymptotics (32j . This is in contrast with the case of the 6j-symbol 
we have been studying here. The lOj-symbol is thus a good example to illustrate the difference between the reasoning 
of the Section [ill B ll in the large spin limit and the exact treatment (when available). In particular, we do not expect 
the lOj-symbol to satisfy the equation (PT|) at the leading order in the asymptotics. 



IV. TOWARDS RECURSION RELATIONS FOR SPINFOAM AMPLITUDES 



The previous sections were dedicated to the new recursion from the point of view of the 6j-symbol itself. Here we 
take a different perspective and try to sketch some implications for the Ponzano- Regge model [g, 0| . 

An important difference between the standard recursion (1151) and the new one ([7]) that we want to emphasize is 
that the latter acts on the square of the 6j-symbol, as opposed to the 6j-symbol itself. While the 6j-symbol is the 
building block of the Ponzano-Regge model, which is attached to each tetrahedron of a triangulation, its square is 
the amplitude for the 3-sphere triangulated by two tetrahedra which are identified along their boundary, 

^pk(^^) -lt[d9aU ^i9a9b') = E [n(2.- + 1)] {£ £ ■ (26) 

This equation is only formal since it is actually divergent (there is one redundant Dirac delta in the product over all 
a < b). Making sense of that expression normally involves a gauge-fixing process, described in Q, which amounts to 
setting one spin in the summand to zero. 

Our new recursion allows a different analysis (which we unfortunately cannot make complete). Note that the 
standard recursion on the 6j-symbol (jlSp cannot be used in the above formula, since it requires to transform indepen- 
dently each symbol. Hence it is not a consistent symmetry of the amplitude. The new recursion, however, is directly 
a symmetry of the summand. 

The usual recursion (1151) encodes the gauge invariance of the model. From the spacetime point of view, that comes 
from specializing the Biedenharn-Elliott identity which states the invariance under the 2-3 Pachner move. In the 
Hamiltonian framework, it is a quantization of a generator of a gauge symmetry which enables to solve for the wave- 
function on the boundary of a tetrahedron [25j . So it may be tempting to interpret the new recursion as being also 
related to that gauge symmetry. However, that is unclear to us, and we have looked at the algebra of the recursions 
acting on different spins without finding pieces of evidence. 

Instead, the new recursion can be interpreted as a global symmetry of the summand. By global, we mean that 
it acts on all spins, none of them being left unchanged. Together with the result of the previous section, that gives 



^ The lOj-symbol has only ten quantum numbers which do not always determine a 4-simplex. But it is known that there is a regime 
where it does [3l| . 
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the following picture. The shifts on the squared 6j-synibol correspond to a deformation of the triangulation of the 
3-sphere equipped with a discrete metric (given by the arguments of the symbol). Invariance under that deformation 
implies that the metric has to close the boundaries of the tetrahedra (which are identified). Thus, it provides with a 
consistent metric deformation of the triangulation. 

We would like to propose to develop this perspective and write recursion relations for more general triangulations. 
Let us start with a triangulation A of the 3-sphere. We consider a connected graph F on A, and fix all the spins je on 
the edges e e F to fixed values. Looking at the Ponzano-Regge amplitude for that triangulation as a function of the 
spins (je)eer, we would like to derive a recursion relation under shifts of the spins (je)eGr- To approach the problem, 
it has been shown in 33] that we can write this amplitude as a group integral. 

We start by choosing a framing for the graph F, i.e. we embed it faithfully in a closed surface in A. More precisely 
we choose a set of triangles forming a closed surface S, such that the edges of F are all in S. In general, F is not a 
triangulation of the whole surface E, but still defines a cellular decomposition of E whose faces are delimited by the 
edges of F. We require that the faces are all homeomorphic to the 2-disk. Then the result proved in [33J is that for a 
trivial topologjQ ^ , the amplitude has a simple expression, with group elements (gf) attached to the faces of E, 

^Abeer] - / [n%] n^>(5s(e)5t7e)) ' (28) 

f 

Here s(e) and t{e) denote the two faces of E sharing the edge e. Since characters are invariant under taking the 
inverse of the argument, it does not matter which face we choose as the source or the target for each edge. 

Written as such, this amplitude is the generalization of the integral for the squared 6j-symbol ([5]). If we look at 
the graph dual to F in E, we see that the group elements (gf) are actually attached to the vertices of the dual graph. 
Thus it is the square of the evaluation of the spin network state dual to the graph F on the surface E. 

The next step is to re-express the measure using the class angles 9e of gs{e)gt(^) as basic variables, 

ZAbeer] = j[deef\ fli^e) J] (^^'-^ + 1) , (29) 



where fi is the Jacobian of the change of variables. Its analysis has been started in 3J] where it has been proven 
that fi involves a product of Gram determinants similarly to the case of the squared 6j-symbol. However, the squared 
6j-symbol is a special case where all the angular variables 6e are independent. Generically that is not true because 
the elements gf are unit vectors in M^, so that only four of them can be linearly independent. Typically, the integral 
for the lOj-symbol contains five copies of SU(2) which implies that one among the ten angles is determined by 
the others. That is the reason which leads /i to be a constraint in ([24)1 . 

In the absence of constraints, one can use the method developed in the present paper to get recursion formulae. 
In the presence of a constrained measure /i, it is in principle possible to also apply the method which was used on 
the lOj-symbol, provided that the constraints in /i are polynomials in (cos^e). However, the geometric content of fi 
seems to depend on the connectivity matrix of F, so that at the end of the day, such recursion formulae lack a global, 
unified geometric interpretation. 



V. OUTLOOK 

In this short paper, we have adapted a method outlined in ^23] to derive recursion relations for the square of the 
6j-symbol using its expression as a multiple integral over SU(2). This method gives a very easy way to obtain recursion 
relations, which has a simple algebraic origin in the expression of the integral measure. Moreover, we have shown how 
it is simply related at leading order in the semi-classical limit to the closure of the tetrahedron whose quantum edge 
lengths are given by the spins labeling the 6j-symbol. Finally, we have explained how our method can be generalized 
to derive recursion relation for arbitrary spin foam amplitudes in the Ponzano-Regge model. We believe that this 
opens the door to a few issues and possibilities, which we hope to study in future work: 



* Reference [33[ also deals with the more complicated case of a non-trivial surface topology. With a non-trivial topology defined by the 
genus g of the surface S, we need to associate group elements Uc € SU(2) to each cycle C of the surface. Assuming that each edge 
e G r belongs to a single cycle C(e) (or to none), we have 

ZAljeer] = f O'^S/l / lU'^^Uc] U^,A9s(e)Uc'l^f9^l^) , (27) 

f C e 

where the sign ec{e) records the relative orientation of the edge e with respect to the orientation of the cycle. 
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• Since the leading order of the recursion relation has a simple geometrical interpretation as the closure of the 
corresponding classical tetrahedron, we can hope that the next-to-leading order of our new recursion relation 
have a similar straightforward interpretation. It may also be that its analysis is technically less involved that 
when using the standard recursion as done in (l6j . 

• It would be interesting to pursue this approach and further study the measure of the group integral for arbitrary 
spin foam amplitudes, in order to understand what geometric properties can emerge. 

• The existence of recursion relations in spin foam amplitudes is symptomatic of the existence of a symmetry. The 
standard recursion relation is usually related to the topological invariance of the Ponzano-Regge model under 
Pachner moves. Since our recursion relation can be interpreted as deformations of the space-time triangulation, 
it would be interesting to investigate its relation to discrete diffeomorphisms. 

• The specific invariance of 3d gravity and more generally BF theory in higher dimensions has been recently 
revisited [l^,!!!,!!!]- In these works, one considers deformation properties of the full amplitude, using the notion 
of flat discrete connections, but not in the spin representation. That would be highly valuable to understand 
how the different objects used in these two different approaches actually probe the same geometric features. 
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